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HIGHLIGHTS 


►  An  extensive  set  of  simulation  results  is  presented  for  several  ethylene/H2/02/N2  flames. 

►  Different  turbulent  combustion  models  are  studied  (EDC,  pdf-transport  and  flamelets). 

►  To  our  knowledge,  this  paper  reports  the  first  CFD  modeling  results  of  the  JHC  flames. 

►  Flame  temperature  and  species  are  predicted  using  two  kinetic  mechanisms  (GRI3.0  and  P0L1M1). 
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This  paper  reports  the  results  of  a  numerical  investigation  of  several  turbulent  nonpremixed  ethylene  jet 
flames,  either  undiluted  or  diluted  with  hydrogen,  air  or  nitrogen,  which  have  been  experimentally 
studied  using  a  jet  in  hot  coflow  (JHC)  burner.  The  fuel  jet  issues  into  a  hot  and  highly  diluted  coflow 
at  two  02  levels  (3%  or  9%)  and  a  temperature  of  1100  K.  These  conditions  emulate  those  of  moderate 
or  intense  low  oxygen  dilution  (MILD)  combustion.  The  attention  is  mainly  focused  on  assessing  the 
performance  of  different  models  for  the  turbulence— chemistry  interaction.  The  model  predictions 
are  compared  to  experimental  measurements.  The  Eddy  dissipation  concept  (EDC),  flamelet  and 
PDF  transport  models  were  used  in  combination  with  a  two-equation  turbulence  model  (k— e)  and 
two  different  kinetic  mechanisms  (GRI-Mech  3.0  and  POLIMI).  The  effect  of  the  oxygen  concentration  in 
the  coflow  is  well  captured  by  all  models.  A  modified  version  of  the  EDC  model,  recently  proposed  for  the 
modeling  of  methane  MILD  combustion  conditions,  is  also  used.  This  model  significantly  improves  the 
predictions  of  the  EDC  in  the  case  of  the  ethylene  flames  studied  in  this  paper.  The  agreement  between 
measurements  and  predictions  of  the  flamelet  model  is  generally  not  as  good  as  the  PDF  transport  and 
modified  EDC  models,  whose  predictions  are  in  good  agreement  with  the  measurements,  and  are 
improved  especially  for  the  apparent  liftoff  heights  and  the  peak  flame  temperature.  The  effect  of  the  fuel 
composition  on  the  CH20  formation  is  also  discussed,  showing  generally  good  agreement. 

©  2013  Elsevier  Ltd.  All  rights  reserved. 


1.  Introduction 

Increasing  energy  demand  and  the  concerns  surrounding 
combustion  generated  pollution  requires  new  combustion  tech¬ 
nologies  to  be  sought  that  are  more  efficient  and  environmentally 
sustainable.  MILD  (moderate  or  intense  low  oxygen  dilution) 
combustion  [1]  is  an  innovative  technique  which  offers  a  possible 
solution  to  this  problem  as  it  allows  a  simultaneous  enhancement 
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in  thermal  efficiency  of  combustion  systems  and  lower  pollutant 
emissions.  This  combustion  regime  proceeds  in  an  environment  of 
low  oxygen  concentration  at  high  oxidizer  temperatures  to  create 
a  more  uniform  and  distributed  combustion  region  which  gives 
higher  efficiency,  radiation  flux  and  lower  levels  of  pollutants.  As 
a  result,  the  combustion  process  takes  place  in  a  more  diluted 
fashion  and  thus  the  temperature  peaks  that  are  largely  responsible 
for  NO*  formation  can  be  avoided.  In  comparison  with  conventional 
combustion  regimes,  MILD  combustion  has  been  reported  to 
provide  energy  savings  of  approximately  30%,  a  30%  reduction  in 
furnace  size  and  a  25%  reduction  in  pollutant  emissions  [2]. 

W tinning  and  Wiinning  [3]  showed  MILD  combustion  (also 
referred  to  as  flameless  oxidation,  FLOX®)  to  be  highly  efficient  with 
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low  NO*  emissions.  There  have  been  numerous  MILD  combustion 
studies  in  semi-industrial  scale  furnaces  [4,5]  and  burners  [6—8], 
both  experimentally  and  using  mathematical  modeling  based  on 
the  Eddy  dissipation  concept  (EDC)  [9,10]  and  conditional  moment 
closure  (CMC)  [11  ].  Danon  and  coworkers  recently  studied  a  multi¬ 
burner  flameless  combustion  furnace.  They  first  studied  exper¬ 
imentally  [12]  the  effect  of  the  position  of  the  burners  and  found 
the  optimal  operating  conditions.  Then,  they  performed  numerical 
simulations  using  the  EDC  model  to  achieve  more  insight  regarding 
the  effects  of  burner  positioning  in  their  multi-burner  flameless 
combustion  furnace  [13]. 

A  review  of  the  research  efforts  in  modeling  MILD  combustion 
have  been  presented  by  Aminian  et  al.  [14]  and  De  et  al.  [15].  These 
authors  also  discussed  the  applicability  of  the  Eddy  dissipation 
concept  to  the  modeling  of  jet  in  hot  coflow  (JHC)  flames  and 
proposed  a  modified  EDC  model  capable  of  enhancing  the  pre¬ 
dictions  of  mean  temperature  field  in  JHC  conditions.  The  interest 
in  the  JHC  burner  is  motivated  by  its  ability  to  emulate  MILD 
combustion  under  controlled  conditions. 

Previous  studies  have  mainly  concentrated  on  the  combustion  of 
CH4  or  CH4/H2  mixtures,  while  in  this  study  the  focus  is  on  the 
modeling  of  ethylene  (C2H4)  flames.  The  interest  toward  ethylene  is 
motivated  by  the  significance  of  this  species  as  an  important  inter¬ 
mediate  in  the  oxidation  of  larger  hydrocarbons,  and  also  to  bridge  the 
gap  between  methane  and  the  more  complex  chemistry  of  practical 
fuels.  In  order  to  advance  fundamental  understanding  of  MILD  com¬ 
bustion,  Medwell  et  al.  [16]  studied  the  combustion  of  C2H4/H2/O2/N2 
mixtures  in  the  JHC  burner  of  the  University  of  Adelaide.  In  the  current 
work  we  discuss  the  results  of  CFD  modeling  of  the  JHC  experiments, 
with  particular  attention  to  the  effect  of  the  models  adopted  with 
a  RANS  approach  for  the  modeling  of  turbulence  and  turbulence- 
chemistry  interactions.  In  particular,  we  examine  the  effects  of  com¬ 
bustion  and  turbulence  models,  differential  diffusion  and  chemical 
kinetic  mechanisms  on  the  accuracy  of  the  predictions. 

2.  Burner  description  and  computational  domain 

The  burner  in  this  study  is  the  jet  in  hot  coflow  (JHC)  burner, 
used  previously  by  Medwell  et  al.  [16,17]  and  is  shown  in  Fig.  1.  It 
consists  of  a  central  insulated  fuel  jet  (0  4.6  mm)  within  an  annular 
coflow  (0  82  mm)  of  hot  exhaust  products  from  a  premixed  sec¬ 
ondary  burner  mounted  upstream  of  the  jet  exit  plane.  The  02  level 
of  the  coflow  is  controlled  by  the  secondary  porous  burner  to  give 
O2  levels  of  3  or  9%  (volumetric),  while  the  coflow  temperature  and 
exit  velocity  were  kept  constant  at  1100  K  and  2.3  m/s.  The  fuel 
used  in  the  jet  is  ethylene  (>99%  C2H4),  either  undiluted,  or  diluted 
with  hydrogen  (H2),  air,  or  nitrogen  (N2).  Table  1  shows  the  com¬ 
position  and  operating  conditions  of  the  fuel  and  oxidant  streams 
for  the  different  cases  studied  in  this  paper.  A  description  of  the 
experimental  details  is  provided  in  Ref.  [16]  and  also  in  Ref.  [17] 
which  refers  to  CH4/H2  flames  and  contains  additional  details  on 
the  experimental  techniques  used  for  the  simultaneous  imaging  of 
OH,  CH20,  and  temperature  in  the  same  burner. 

It  is  worth  noting  that  although  the  experimental  results  pre¬ 
sented  in  the  present  paper  are  shown  on  the  graphs  only  as  a  line, 
there  are  inherent  experimental  errors  and  uncertainties.  Factors 
such  as  correction  for  the  Rayleigh  cross-section  (for  temperature), 
collisional  quenching  (for  OH),  laser  power  fluctuations,  and  noise 
all  affect  the  experimental  measurements.  Analysis  of  measure¬ 
ments  in  a  flat-flame  indicates  systematic  errors  of  5%  for  the  OH 
and  2.5%  for  temperature.  In  a  turbulent  flame  these  errors  are 
estimated  to  be  less  than  20%  in  the  worst  case,  and  for  the  aver¬ 
aged  data  presented  in  this  paper  less  than  10%. 

In  this  work,  numerical  simulations  are  performed  to  model  the 
flames  in  the  JHC  burner  developed  and  experimentally  studied  by 


Medwell  et  al.  [16].  The  flames  were  simulated  with  the  commer¬ 
cial  CFD  code  FLUENT  12.1.4.  A  two-dimensional  steady-state  sim¬ 
ulation  of  the  physical  domain  was  considered  due  to  the  symmetry 
of  the  system.  Fig.  2  shows  the  computational  grid  used  to  simulate 
the  flames  and  the  boundary  conditions  adopted.  The  grid  domain 
is  400  mm  in  the  axial  direction  and  120  mm  in  the  radial  direction 
from  the  jet  exit.  It  is  a  structured  non-uniform  grid  with  34,200 
cells,  which  provides  high  resolution  in  the  reaction  zone  region 
and  close  to  the  inlets  and  save  computational  effort  elsewhere. 
This  grid  layout  was  selected  based  on  an  examination  of  different 
cell  sizes,  and  no  further  significant  change  was  obtained  for  finer 
cells.  The  flame  is  non-confined  and  hence  a  pressure  outlet 
condition  is  used  (see  Fig.  2)  assuming  ambient  air  backflow  con¬ 
ditions.  The  simulations  were  performed  for  two  different  O2  mole 
fractions  in  the  coflow  stream,  namely  3%  and  9%.  Fig.  3  anticipates 
the  results  of  some  model  predictions  which  are  discussed  later.  It 
compares  predicted  temperature  fields  and  photographs  of  the 
flame  [16]  and  clearly  shows  the  effect  of  the  different  amount  of 
oxygen  in  the  oxidizer.  Consistent  with  the  ‘flameless’  character¬ 
istics  of  MILD  combustion,  the  photographs  may  indicate  that  the 
flames  appear  lifted,  but  reactions  are  indeed  present  between  the 
jet  exit  plane  and  the  visible  flame.  The  presence  of  soot  at  the 
downstream  locations  is  a  result  of  the  entrainment  of  surrounding 
air  beyond  the  confines  of  the  coflow,  where  MILD  combustion 
conditions  are  no  longer  maintained.  As  expected,  combustion 
occurs  with  weaker  or  higher  intensity  according  to  the  oxygen 
content  in  the  oxidizer  stream. 

The  aim  of  this  work  is  first  to  identify  key  modeling  issues  for 
predicting  temperature  and  chemical  species  in  the  ethylene  JHC 
experiments  and  then  to  discuss  the  effect  of  fuel  composition 
on  the  reactivity  of  the  reacting  mixture  using  two  kinetic 
mechanisms. 


540 


S.R.  Shabanian  et  al.  /  Applied  Thermal  Engineering  52  (2013)  538-554 


Table  1 

Operating  conditions  of  two  inlet  streams  (on  a  molar  basis)  for  the  different  case  studies. 


Flame 

Fuel  jet 

Oxidant  flow 

U  (m/s) 

T(I<) 

c2h4 

H  2 

o2 

N2 

U  (m/s) 

T(K) 

02 

n 

O 

h2o 

n2 

C2H4— 3%  02 

17.5 

305 

100 

0 

0 

0 

2.3 

1100 

3 

3 

10 

84 

C2H4— 9%  02 

17.5 

305 

100 

0 

0 

0 

2.3 

1100 

9 

3 

10 

78 

C2H4/H2— 3%  o2 

30.6 

305 

50 

50 

0 

0 

2.3 

1100 

3 

3 

10 

84 

C2H4/H2— 9%  02 

30.6 

305 

50 

50 

0 

0 

2.3 

1100 

9 

3 

10 

78 

C2H4/air— 3%  02 

27.3 

305 

25 

0 

15.75 

59.25 

2.3 

1100 

3 

3 

10 

84 

C2H4/air— 9%  02 

27.3 

305 

25 

0 

15.75 

59.25 

2.3 

1100 

9 

3 

10 

78 

C2H4/N2— 3%  02 

27.3 

305 

25 

0 

0 

75 

2.3 

1100 

3 

3 

10 

84 

C2H4/N2— 9%  02 

27.3 

305 

25 

0 

0 

75 

2.3 

1100 

9 

3 

10 

78 

3.  Modeling 

Christo  and  Dally  [9]  have  shown  that  in  the  case  of  the  JHC 
burner  fed  with  CH4/H2,  differential  diffusion  effects  have  a  strong 
influence  on  the  accuracy  of  the  predictions  and  that  conserved 
scalar  based  models,  i.e.  the  £/PDF  and  flamelet  models  (adopting 
equilibrium  or  finite  rate  chemistry),  are  inadequate  for  modeling 
jet  in  hot  coflow  (JHC)  flames.  Moreover,  they  also  discussed  the 
importance  of  using  detailed  chemical  kinetics,  rather  than  global 
or  simplified  mechanisms. 

It  is  important  to  notice  that,  compared  to  traditional  diffusion 
flames,  the  uniform  conditions  and  low  oxygen  concentration  in 
MILD  combustion  regime,  lead  to  slower  reaction  rates  and  enhance 
the  influence  of  molecular  diffusion  on  flame  characteristics.  These 
two  effects  challenge  the  applicability  of  simple  combustion  models 
that  assume  fast  chemistry  and  neglect  the  effects  of  differential 
diffusion  [9], 

The  objective  of  the  study  reported  in  this  work  is  to  assess  the 
performance  of  the  different  turbulence— chemistry  interaction 
models  in  combination  with  different  turbulence  models  and 
chemical  kinetic  mechanisms  in  predicting  JHC  ethylene  flames  by 
comparing  predictions  with  experimental  data.  One  of  the  main 
differences,  compared  to  the  previous  methane  (and  CH4/H2)  cases, 
is  that  ethylene  is  characterized  by  a  significantly  higher  reactivity 
and  relatively  more  complex  chemistry.  The  kinetic  model  adopted 
should  be  able  to  properly  characterize  the  role  of  vinyl  radicals  and 
the  higher  acetylene  formation  in  ethylene  combustion.  Moreover, 
ethylene  has  a  greater  tendency  to  form  soot. 

3.1.  Turbulence  and  radiation  models 

The  burner  was  modeled  using  several  turbulence  and  radiation 
models  to  simulate  the  flames.  To  model  turbulence,  three  different 
models  have  been  evaluated:  the  modified  standard  k—e  (modified 


SKE),  the  renormalization  group  (RNG)  k—e  and  the  Reynolds  stress 
model  (RSM).  The  standard  k—e  model  is  a  semi-empirical  approx¬ 
imation  based  on  solving  transport  equations  for  the  turbulence 
kinetic  energy  (k)  and  its  dissipation  rate  (e),  with  both  equations 
containing  a  number  of  empirical  constants.  In  this  model,  the  flow  is 
assumed  fully  turbulent,  and  the  effects  of  molecular  viscosity 
are  neglected.  Assumption  of  an  isotropic  scalar  quantity  for  the  tur¬ 
bulent  viscosity  is  another  approximation  of  the  standard  k—e  model. 
This  model  is  known  for  its  shortcomings  in  predicting  round  jets.  In 
particular,  the  standard  k—e  model  over-predicts  the  decay  rate  and 
the  spreading  rate  of  a  round  jet.  The  Cit  parameter  is  modified  to 
a  value  of  1.60  for  self-similar  round  jets  [  18].  The  RNG  k—e  model  has  a 
formulation  similar  to  that  of  the  SKE  model,  although  it  uses  different 
constants  and  has  additional  terms  and  functions  in  the  transport 
equations  of  k  and  e.  In  this  model,  there  are  exact  expressions  for  the 
turbulent  Prandtl  and  Schmidt  numbers.  The  RNG  model  is  expected 
to  yield  better  results  in  predicting  round  jets,  swirling  flows  and 
flows  with  a  recirculation  region  in  comparison  with  the  SKE  model 
[19].  The  Reynolds  stress  model  (RSM)  involves  calculation  of  the 
individual  Reynolds  stresses  using  differential  transport  equations. 
The  Navier— Stokes  equations  with  Reynolds  stress  model  become 
stiff  and  time  consuming  in  convergence  compared  to  the  k—e  model. 
The  RSM  does  not  use  the  Boussinesq  hypothesis  and  isotropic 
turbulent  viscosity  assumption.  Gran  et  al.  [20]  have  reported 
improvement  in  prediction  of  observed  flow  patterns  in  the  near 
burner  zone  for  combusting  flows  using  the  RSM. 

As  suggested  by  Christo  and  Dally  [9],  Frassoldati  et  al.  [10],  and 
Aminian  et  al.  [14],  the  discrete  ordinate  (DO)  radiation  model  [21] 
has  been  used  to  calculate  the  flame  thermal  radiation.  This  model 
can  be  used  for  a  wide  range  of  optical  thicknesses.  The  model  solves 
a  radiative  transfer  equation  for  a  number  of  discrete  solid  angles 
across  the  computational  domain.  For  the  model,  theta-divisions  and 
phi-divisions  define  the  number  of  control  angles  used  to  discretize 
each  octant  of  the  angular  space.  Also,  theta-pixels  and  phi-pixels  are 
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Coflow — ► 

Axis 


Fig.  2.  Two-dimensional  view  of  the  computational  domain  and  boundary  conditions.  Boundary  conditions  are  indicated  in  the  figure.  Walls  are  assumed  adiabatic. 
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Fig.  3.  Experimental  photographs  [16]  and  predicted  temperature  field  (the  modified  EDC  combustion  model,  modified  SKE  turbulence  model,  GRI-Mech  3.0  kinetic  mechanism  are 
used  in  generating  the  temperature  field)  of  pure  ethylene  flames  at  i?ejet  =  10,000  for  two  coflow  02  levels.  Dotted  lines  show  the  locations  of  the  measurements,  which  were  taken 
along  the  axial  height  35  and  125  mm  from  the  inlet. 


used  to  control  the  pixelation  that  accounts  for  any  control  volume 
overhang.  The  model  computes  the  absorption  coefficient  from  the 
weighted  sum  of  the  gray  gas  model  (WSGGM)  in  which  spatial 
variation  in  the  total  emissivity  is  computed  as  a  function  of  gas 
composition  and  temperature.  The  WSGGM  is  a  model  with 
reasonable  accuracy  that  takes  into  account  particular  absorption 
bands.  It  is  worth  noting  that  the  effects  of  radiation  on  the  flame  are 
negligible  when  in  the  confines  of  the  hot  coflow  (axial  distance 
below  ~120  mm)  (refer  to  Fig.  SI  in  the  supplemental  material  for 
further  details). 


3.2.  Turbulent  combustion  models 

The  main  challenges  in  modeling  turbulent  combustion  are 
handling  the  mean  rate  of  reaction  and  adequate  representation  of 
the  chemistry  in  the  model.  In  this  study  the  capability  of  three 
classes  of  combustion  models  for  predicting  MILD  flames  is  assessed: 
volumetric  reaction-based  models,  conserved  scalar  based  models 
and  PDF  transport  models.  The  Eddy  dissipation  concept  (EDC) 
model  [22,23]  and  the  modified  EDC  [15,24]  belong  to  the  first 
category,  the  Steady  Laminar  Flamelets  (SLF)  model  [25,26]  to  the 
second.  Finally,  the  composition  PDF  transport  models  allow  simu¬ 
lating  finite-rate  chemistry  in  turbulent  flames  without  the 
assumptions  needed  to  model  the  unknown  terms  which  arise  from 
the  Reynolds-averaging  of  the  energy  and  species  balance  equations. 
Nevertheless,  the  PDF  transport  models  require  the  modeling  of  the 
molecular  mixing  of  species  and  heat  which  is  usually  the  most 
significant  source  of  error  in  the  PDF  transport  approach. 

3.2.3.  Eddy  dissipation  concept  model 

The  Eddy  dissipation  concept  (EDC)  model  has  the  advantage  of 
incorporating  detailed  kinetics  at  a  computational  cost  which  is 
quite  moderate  compared  to  more  advanced  models  such  as  the 
transported  PDF  method,  especially  if  the  in-situ  adaptive  tabu¬ 
lation  (ISAT)  model  [27]  is  used.  As  already  discussed  by  De  et  al. 
[15],  however,  this  advantage  comes  at  the  cost  of  a  less  accurate 
description  of  turbulent  fluctuations. 

In  a  turbulent  environment,  combustion  takes  place  where  there 
is  a  molecular  mixing,  i.e.  at  small  turbulence  scales.  According  to 
the  EDC  model,  the  chemical  reactions  occur  only  in  small  scale 
micro-mixed  turbulent  structures  known  as  fine  structures. 


The  evolution  of  species  mass  fractions  inside  the  fine  scales  (w*) 
is  then  computed  by  integrating  the  chemistry  within  these  fine 
scales  for  a  residence  time  z*. 


(1) 


where  CT  is  a  time  scale  constant  equal  to  0.4082  in  the  standard 
EDC  model  and  becomes  Cr  =  3.0  in  the  modified  EDC  proposed  by 
De  et  al.  [15],  Here  v  is  the  laminar  kinematic  viscosity. 

The  size  of  the  fine  scales  is  modeled  as  a  function  of  the  local 
turbulent  properties: 


(2) 


v*  =  n-v 


(3) 


The  constant  Cy  is  equal  to  2.1377  [22]. 

Based  on  the  mass  transfer  between  the  fine  structures  and  their 
surroundings,  the  mean  reaction  term  becomes: 


Ri  = 


pyj 

**(i  -r*3) 


(w*  -  w?) 


(4) 


where  p  is  the  density  and  w,  is  the  mass  fraction  of  each  species. 
The  EDC  model  can  incorporate  detailed  chemical  mechanisms  into 
turbulent  reacting  flows  and  can  be  used  when  the  assumption  of 
fast  chemistry  is  invalid.  However,  typical  mechanisms  are  invar¬ 
iably  stiff  and  their  numerical  integration  is  computationally  costly. 


3.2.2.  Steady  laminar  flamelet  model 

The  basic  assumption  is  that  instantaneous  thermo— chemical 
state  of  the  fluid  is  related  to  a  conserved  scalar  quantity  known 
as  the  mixture  fraction.  In  this  way  the  species  transport  equations 
can  be  reduced  to  a  transport  equation  for  the  mixture  fraction  ? 
and  one  for  its  variance  f'2: 


£<rf)  +  £(,u>f) 


dXj  \Sct  dXj] 


(5) 
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(6) 


The  constant  C0  inside  the  dissipation  term  can  be  derived  by 
turbulent  spectral  analysis  and  is  usually  set  at  2.  If  the  system  is 
not  adiabatic,  the  enthalpy  balance  equation  must  be  solved: 


3 

at 


( PH)  + 


a  fktdH\ 
dXj  \CpdxJ 


+  SH 


(7) 


where  kt  and  Cp  are  the  thermal  conductivity  and  specific  heat  of 
the  mixture,  and  SH  is  a  generic  source  term  which  accounts  for  the 
non-adiabatic  behavior  of  the  system. 

The  temperature  and  thermo— chemical  variables  are  extracted 
from  a  flamelet  library,  in  which  the  temperature  and  composition 
corresponding  to  each  value  of  the  mean  mixture  fraction,  mixture 
fraction  variance  and  enthalpy  are  stored.  The  average  value  of  the 
generic  scalar  \j/  stored  in  the  library  is  evaluated  by  the  following 
integral: 


i=  f  ftOO-KlH)  d?  (8) 

0 


where  /^(f)  is  the  mixture  fraction  PDF  (probability  distribution 
function)  and  represents  the  relationship  between  mixture 

fraction  and  enthalpy  to  the  scalar  \j/.  In  this  work  a  (3-PDF  is 
employed.  The  function  which  takes  into  account  the 

treatment  of  complex  chemistry,  can  be  modeled  following  a  dif¬ 
ferent  approach.  In  the  case  of  the  steady  laminar  flamelet  model 
the  turbulent  flame  is  considered  as  an  ensemble  of  discrete,  steady 
laminar  flames,  called  flamelets. 


3.2.3.  Composition  PDF  transport  model 

The  composition  PDF  transport  model  allows  simulating  finite- 
rate  chemical  kinetic  effects  in  turbulent  reacting  flows  as  well  as 
flame  extinction  and  ignition.  These  PDF  transport  simulations  are 
computationally  expensive,  but  the  computational  power  currently 
available  allows  the  use  of  such  an  approach  for  two-dimensional 
simulations. 

An  alternative  to  Reynolds-averaging  the  species  and  energy 
equations  (which  leads  to  unknown  terms  for  the  turbulent  scalar 
flux  and  the  mean  reaction  rate  that  needs  to  be  modeled)  is  to 
derive  a  transport  equation  for  their  single-point,  joint  probability 
density  function  (PDF)  [28].  The  main  strength  of  the  PDF  transport 
approach  is  that  the  highly-nonlinear  reaction  term  is  completely 
closed  and  does  not  require  modeling.  This  PDF,  denoted  by  P,  can 
be  considered  to  be  proportional  to  the  fraction  of  the  time  that  the 
fluid  spends  at  each  species  and  temperature  state:  P  has  N  +  1 
dimensions  for  the  N  species  and  temperature.  From  the  PDF,  any 
thermo— chemical  moment  (e.g.  mean  reaction  rate  or  mean  or  RMS 
temperature)  can  be  calculated.  The  composition  PDF  transport 
equation  is  derived  from  the  Navier— Stokes  equations  as  [29]: 


at 


<«+A(m.p)  +  _i_(pSlP)  _ 


+ 


a 

Wk  L' 


,  1  9/i,k 

^  p  dXi 


*)P 


(9) 


where  P  is  the  Favre  joint  PDF  of  composition,  p  is  the  mean 
density  of  the  fluid,  u,-  is  the  Favre  mean  fluid  velocity  vector,  S ^  is 
the  reaction  rate  for  species  k,  is  the  composition  space  vector, 
uj’is  the  fluid  velocity  fluctuation  vector,  and  is  the  molecular 


diffusion  flux  vector.  The  generic  notation  <A|B>  indicates  the 
conditional  probability  of  event  A,  given  the  occurrence  of  event  B. 
In  Equation  (1 ),  the  terms  on  the  left-hand  side  are  closed,  whereas 
those  on  the  right-hand  side  are  not  and  require  modeling.  The 
turbulent  scalar  flux  term  is  modeled  by  the  gradient  diffusion 
assumption: 


where  pt  is  the  turbulent  viscosity  and  Sct  is  the  turbulent  Schmidt 
number. 

In  PDF  methods,  the  effect  of  molecular  diffusion  on  the  compo¬ 
sition  is  represented  by  a  mixing  model.  Two  mixing  models,  the 
modified  curl  (MC)  model  [30,31]  and  Euclidean  minimum  spanning 
tree  (EMST)  model  [32]  were  used  in  this  study  to  model  mixing  at 
the  smallest  molecular  scales.  For  the  MC  model,  a  few  particle  pairs 
are  selected  at  random  from  all  the  particles  in  a  cell,  and  their 
individual  compositions  are  moved  toward  their  mean  composition. 
The  number  of  particle  pairs  selected  is  calculated  from  the  algorithm 
according  to  the  total  number  of  particles  and  the  local  turbulent  time 
scale  (fc/e).  For  each  particle  pair,  a  uniform  random  number  /?  is 
selected  and  each  particle’s  composition  0  is  moved  toward  the  pair’s 
mean  composition  by  a  factor  proportional  to  fk 


0}  =  (1  -  8)0?  +  8 


(0°ml  +  0?mj 
(m,-  +  rrij) 


0, 


.  n  K0?mi  +  0?mj 

5  =  (1  -  0)0?  +  8^-. - 

]  1  ( mj  +  ntj ) 


(H) 


(12) 


where  0;  and  0j  are  the  composition  vectors  of  particles  i  and  j,  and 
m;  and  m,  are  the  masses  of  particles  i  and  j. 

Physically,  mixing  occurs  between  fluid  particles  that  are 
adjacent  to  each  other.  The  MC  mixing  model  takes  no  account  of 
this  localness,  which  can  be  a  source  of  error.  The  advantage  of  the 
Euclidean  minimum  spanning  tree  (EMST)  model  is  that  it  mixes 
particle  pairs  that  are  close  to  each  other  in  composition  space.  In 
this  way  the  model  takes  into  account  the  physical  constraint  that 
mixing  occurs  between  fluid  particles  that  are  adjacent  to  each 
other.  Since  scalar  fields  are  locally  smooth,  particles  that  are  close 
in  composition  space  are  likely  to  be  close  in  physical  space.  The 
particle  pairing  is  determined  by  a  Euclidean  minimum  spanning 
tree,  which  is  the  minimum  length  of  the  set  of  edges  connecting 
one  particle  to  at  least  one  other  particle.  The  EMST  mixing  model  is 
more  accurate  than  the  MC  and  interaction  by  exchange  with  the 
mean  (IEM)  mixing  models,  but  is  slightly  more  computationally 
expensive.  Details  on  the  EMST  model  can  be  found  in  Ref.  [32]. 


3.3.  Kinetic  mechanisms 

Due  to  the  recommendations  of  previous  studies  [9,10,14],  here 
we  adopt  only  detailed  kinetic  mechanisms.  Two  detailed  schemes, 
the  GRI-Mech  3.0  [33]  and  POLIMI  ethylene  mechanism  [34]  were 
evaluated.  The  GRI-Mech  3.0  mechanism  was  employed  without 
the  NO*  reactions,  resulting  in  219  reversible  reactions  relating  36 
chemical  species.  The  POLIMI  ethylene  mechanism  used  in  this 
work  is  a  reduced  version  of  the  combustion  mechanism  developed 
at  Politecnico  di  Milano.  The  complete  scheme  is  available  on  the 
website  indicated  in  ref  [34].  The  kinetic  model  used  in  this  study 
was  described  and  validated  before  using  a  large  set  of  exper¬ 
imental  data  [35—40],  including  MILD  combustion  conditions 
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[10,41].  We  recently  further  validated  the  kinetic  mechanism  for 
ethylene  using  the  coflow  flame  measurements  of  McEnally  and 
Pfefferle  [42],  as  discussed  by  Cuoci  et  al.  [43],  The  numerical  results 
were  found  to  be  in  satisfactory  agreement  with  the  experimental 
data,  demonstrating  the  reliability  of  the  kinetic  mechanism.  Due  to 
the  large  dimensions  of  this  scheme,  it  is  not  suitable  for  CFD 
simulations  of  turbulent  flames.  Therefore,  a  reduced  version  was 
obtained  using  an  algorithm  [44]  which  combines  the  reaction  path 
analysis  and  the  directed  relation  graph  method  [45].  The  reduced 
model  contains  24  species  involved  in  155  reactions  and  is  able  to 
reproduce  the  performance  of  the  original  mechanism  over  a  wide 
range  of  operating  conditions  (0.5  <  4>  <  2.0)  at  high  temperatures 
(above  1000  I<)  and  atmospheric  pressure. 

4.  Comparison  with  experimental  measurements 

4.2.  Effect  of  turbulence  model 

Turbulent  mixing  has  a  significant  effect  on  the  flow  field, 
temperature  and  turbulence— chemistry  interaction.  In  this  study, 
turbulence  was  modeled  via  the  RANS  approach,  using  the  modified 
standard  k—e  (modified  SKE),  RNG  k—e  and  the  Reynolds  Stress 
Model  (RSM).  Using  the  modified  EDC  combustion  model  [15], 
discrete  ordinates  (DO)  radiation  model  and  GRI-Mech  3.0  mech¬ 
anism,  calculations  for  the  3%  and  9%  02  flames  were  performed. 
For  comparing  the  different  turbulence  models  in  this  section, 
only  the  GRI-Mech  3.0  mechanism  was  used.  Fig.  4  shows  radial 
distributions  of  the  temperature  obtained  using  the  three  turbu¬ 
lence  models  at  an  axial  location  of  35  mm  above  the  jet  exit  plane 
for  the  3%  02  flames.  It  can  be  observed  that  discrepancies  between 
the  experiments  and  model  predictions  are  larger  for  the  RNG  k—e 
model.  It  is  noteworthy  that  the  RNG  model  fails  to  predict  a  reac¬ 
tion  for  all  but  the  pure  ethylene  flame.  The  poor  agreement  of  the 


C2H4-3%  02-axial  location  35  mm 


C2H4/H2-3%  02-axial  location  35  mm 


RNG  model  was  also  observed  by  Christo  and  Dally  [9],  Frassoldati 
et  al.  [10]  and  De  et  al.  [15]  for  methane  and  methane/hydrogen  jet 
flames.  Moreover,  the  predictions  of  the  modified  SKE  and  RSM 
models  are  almost  similar.  This  figure  shows  generally  good 
agreement,  and  excellent  agreement  between  predicted  profiles  by 
the  modified  SKE  and  RSM  with  experiments  for  the  C2H4/H2  and 
C2H4/N2  flames. 

The  temperature  predictions  at  the  axial  location  of  125  mm  are 
more  sensitive  to  the  turbulence  model  used.  Nonetheless,  the 
general  agreement  between  the  models  and  the  experiments  for 
the  3%  02  flames  is  good  (refer  to  Fig.  S2  (Supplementary  data)). 
At  this  location  it  is  noted  that  the  measured  high  temperature 
peak  in  the  C2H4/air  flame  is  not  captured  in  the  models,  which 
could  be  a  consequence  of  the  complex  three-stream  mixing  (due  to 
the  entrainment  of  surrounding  room  air).  Importantly,  the  higher 
centerline  temperature  for  the  air-diluted  flame  is  in  agreement 
between  the  experiments  and  models.  For  the  9%  02  flames,  the 
RSM  and  modified  SKE  predict  the  location  of  maximum  temper¬ 
ature  better  than  RNG  k—e  model  (refer  to  Fig.  S3  (Supplementary 
data)).  In  addition,  the  RSM  model  predictions  have  smaller 
deviations  from  the  experiments  in  comparison  to  the  modified  SKE 
predictions.  In  most  cases,  however,  the  measured  peak  tempera¬ 
ture  is  significantly  lower  than  that  from  the  models.  This  effect  is 
believed  to  be  an  artifact  of  the  lifted  characteristics  of  the  exper¬ 
imental  flames  that  has  not  been  captured  by  the  models.  The 
comparison  between  predicted  and  measured  [16]  OH  number 
densities  allows  better  assessment  of  the  effect  of  the  different 
turbulence  models.  These  results  (not  shown)  demonstrated  that 
the  modified  SKE  and  RSM  models  are  more  accurate  than  the  RNG 
k—e  model  in  all  these  flames.  In  this  study,  the  modified  SKE  model 
is  then  used  in  all  of  the  simulations  in  order  to  avoid  the  additional 
computational  expense  of  the  Reynolds  stress  model,  and  to  con¬ 
centrate  on  the  study  of  turbulence— chemistiy  interactions. 
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Fig.  4.  Effect  of  turbulence  models  (modified  EDC  [15]  and  GRI-Mech  3.0  mechanism  used  for  the  combustion  modeling  and  discrete  ordinates  (DO)  model  for  radiation):  com¬ 
parison  between  experimental  measurement  and  model  predictions  at  axial  location  of  35  mm  for  the  3%  02  coflow  oxidant  stream  and  different  fuels:  (a)  C2H4,  (b)  C2ELi/air,  (c) 
C2H4/H2,  (d)  C2H4/N2. 
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4.2.  Effect  of  turbulent  combustion  models 

In  this  section  the  effect  of  the  models  adopted  to  describe  the 
interactions  between  turbulence  and  chemistry  is  discussed. 
Experimental  measurements  [16]  are  compared  to  predictions 
obtained  using  the  models  described  in  Section  3.2.  Five  turbulent 
combustion  models  are  used  to  simulate  the  ethylene  flames  in  the 
JHC  burner  using  the  GRI-Mech  3.0  kinetic  mechanism. 

4.2 A.  Ethylene  flames 

Fig.  5  illustrates  the  predictions  of  the  different  turbulent 
combustion  models  for  the  3%  02  and  9%  02  ethylene  flames,  with 
measured  values  of  temperature  and  OH  number  density  radial 
profiles  at  an  axial  location  of  35  mm  and  125  mm.  It  shows  that  all 
the  models  but  EDC  give  good  temperature  predictions  for  the  3% 
02  flame  at  the  35  mm  location  (Fig.  5a).  Fig.  5b  shows  a  reasonable 
agreement  between  predicted  and  measured  OH  number  density 
radial  profile  when  the  EDC  model  is  used  and  some  discrepancies 
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between  the  other  models  and  experimental  values.  The  higher  OH 
values  of  EDC  model  are  expected  because  of  the  higher  predicted 
value  of  temperature  at  the  35  mm  location.  It  should  be  recognized 
that  there  is  some  uncertainty  in  the  absolute  values  of  the  meas¬ 
ured  OH  number  density,  but  importantly  the  location  and  features 
of  the  OH  profiles  are  consistent  between  the  experiments  and  the 
models.  Fig.  5b  also  shows  that  the  PDF  transport  models  (modified 
curl  and  EMST)  predict  noticeable  equilibrium  OH  concentration  in 
the  coflow  region,  while  the  other  models  predict  almost  zero  value 
for  OH  in  the  coflow  stream.  At  the  125  mm  axial  location,  as 
observed  in  Fig.  5c,  the  models  of  modified  EDC,  flamelets,  PDF 
transport  modified  curl  and  EMST  predict  the  temperature  profile 
reasonably  well,  but  the  EDC  model  noticeably  over-predicts  the 
temperature  peak.  In  addition,  noticeable  discrepancies  with 
measured  values  of  OH  number  density  can  be  observed  for  all  the 
used  models  as  shown  in  Fig.  5d.  This  could  be  a  result  of  the 
complex  three-stream  mixing  with  the  surrounding  air  that  is 
difficult  to  incorporate  into  the  models,  as  already  noted  by  Christo 
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Fig.  5.  Effect  of  turbulent  combustion  models  (adopting  the  modified  SKE  turbulence  model,  GRI-Mech  3.0  mechanism  and  (DO)  radiation  model):  comparison  between  exper¬ 
imental  measurement  [16]  and  model  predictions  (temperature  (a,  c,  e)  and  OH  number  density  (b,  d,  f))  for  3%  02  and  9%  02  ethylene  flames. 
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and  Dally  [9],  For  the  9%  02  flames,  all  models  tend  to  overestimate 
the  magnitude  of  temperature  peak.  This  discrepancy  is  likely 
a  result  of  the  lifted  flame  behavior  that  was  observed  in  the 
experiments  which  has  not  been  replicated  in  the  models.  Inter¬ 
estingly,  all  models  generally  provide  good  agreement  with 
experiment.  Fig.  5f  reveals  that  the  EDC  model  can  predict  accu¬ 
rately  the  magnitude  of  OF!  number  density  while  the  flamelets  and 
PDF  transport  models  give  better  predictions  of  the  peak  location  of 
OH  concentration. 

4.2.2.  Ethylene/hydrogen  flames 

A  comparison  between  the  modeling  and  experimental  data 
is  presented  in  Fig.  6  for  the  two  ethylene— hydrogen  flames  (3%  02 
and  9%  02  in  the  coflow)  and  two  axial  locations.  The  figure  shows 
that  for  the  3%  02  flame,  the  predicted  temperature  of  the  modified 
EDC,  PDF  transport  modified  curl  and  EMST  models  at  35  mm  and 
125  mm  are  in  good  agreement  with  experimental  measurements. 
Deviations  in  the  absolute  values  of  OH  profiles  from  the  experi¬ 
ments  are  quite  large  for  the  all  models  used,  though  the  locations 


of  the  peaks  are  in  agreement.  For  the  9%  02  flame,  the  results  of  all 
models  except  EDC  show  reasonable  agreement  in  comparison  to 
experimental  temperature  measurements.  In  addition,  the  OH 
number  density  values  in  the  9%  02  flame  have  comparable  mag¬ 
nitude  with  experimental  measurements  for  all  of  the  models. 
Again,  similar  to  the  observation  for  the  pure  ethylene  flames,  the 
conserved  scalar  based  models  give  a  better  prediction  of  the  peak 
location  of  OH  concentration. 

4.2.3.  Ethylene/air  flames 

Fig.  7  shows  the  temperature  and  OH  number  density  profiles  of 
3%  and  9%  02  flames  at  35  mm  and  125  mm  for  the  ethylene— air 
fuel  case.  As  can  be  seen  from  Fig.  7a,  the  modified  EDC  and  PDF 
transport  EMST  give  satisfactory  predictions  of  temperature  at 
35  mm  for  the  3%  02  flame.  In  addition,  the  EDC  model  provides 
good  prediction  of  magnitude  of  temperature  peak,  but  the  shape 
(width)  of  temperature  profile  is  not  fully  consistent  with  exper¬ 
imental  trend.  In  this  partially  premixed  case,  it  can  be  also 
observed  that  the  flamelets  model  is  too  reactive  in  the  region  close 
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Fig.  6.  Effect  of  turbulent  combustion  models  (adopting  the  modified  SKE  turbulence  model,  GRI-Mech  3.0  mechanism  and  (DO)  radiation  model):  comparison  between  exper¬ 
imental  measurement  [16]  and  model  predictions  (temperature  (a,  c,  e)  and  OH  number  density  (b,  d,  f))  for  3%  02  and  9%  02  ethylene/hydrogen  flames. 
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Fig.  7.  Effect  of  turbulent  combustion  models  (adopting  the  modified  SKE  turbulence  model,  GRI-Mech  3.0  mechanism  and  (DO)  radiation  model):  comparison  between  exper¬ 
imental  measurement  [16]  and  model  predictions  (temperature  (a,  c,  e)  and  OH  number  density  (b,  d,  f))  for  3%  02  and  9%  02  ethylene/air  flames. 


to  the  axis  and  the  prediction  of  this  model  is  not  accurate  at 
location  of  35  mm.  Fig.  7b  reveals  that  the  EDC  and  modified  EDC 
models  provide  reasonable  OH  values  in  3%  02  flame  at  35  mm.  Also 
in  this  case,  the  OH  magnitude  in  the  coflow  stream  for  PDF 
transport  models  is  comparable  to  experimental  values,  while  the 
other  models  predict  much  lower  OH  values  in  the  coflow  stream. 
Fig.  7c  shows  that  the  predicted  temperature  profiles  are  in  rea¬ 
sonable  agreement  for  all  combustion  models,  but  the  shape  of 
the  profiles  are  not  correctly  predicted  when  compared  with  the 
experiments.  All  combustion  models  underestimate  the  OH  num¬ 
ber  density  at  125  mm.  It  is  important  to  note  that  Medwell  et  al. 
[16]  showed  that  the  influence  of  the  coflow  persists  ~100  mm 
downstream  of  the  jet  exit  plane  and  that  beyond  this  point  the 
surrounding  air  begins  to  mix  with  the  jet  and  coflow.  For  this 
reason,  these  discrepancies  may  be  related  to  the  effect  of  the 
turbulence  model  and  especially  the  boundary  conditions  adopted 
for  the  coflow  and  surrounding  air  inlets  (see  Fig.  2).  This  discus¬ 
sion,  already  presented  in  detail  by  Frassoldati  et  al.  [10]  and 
Aminian  et  al.  [14],  is  outside  the  scope  of  this  work,  which  focuses 


mainly  on  the  effect  of  turbulent  combustion  models  on  the  CFD 
predictions.  For  the  9%  02  flame,  it  can  be  seen  that  the  EDC, 
modified  EDC  and  flamelet  models  are  too  reactive  and  over¬ 
predict  the  temperature  and  OH  number  density  at  35  mm. 
Against  these  models,  the  PDF  transport  approach  yield  an  accurate 
prediction  of  temperature  and  OH  profiles  for  the  9%  02  flame, 
particularly  the  EMST  model.  It  can  be  observed  that  the  (apparent) 
liftoff  height  is  very  sensitive  to  the  model  adopted. 

4.2.4.  Ethylene/nitrogen  flames 

Fig.  8  illustrates  the  temperature  and  OH  number  density  pro¬ 
files  of  3%  and  9%  02  flames  at  35  mm  and  125  mm  for  ethylene- 
nitrogen  fuel.  As  expected,  due  to  the  inert  addition  there  is 
a  lower  reactivity  in  the  ethylene— nitrogen  flame  with  respect  to 
other  fuels.  The  apparent  liftoff  height  (not  presented  in  this  paper) 
seems  correctly  predicted  by  all  models.  For  the  3%  02  flame,  all 
combustion  models  give  an  accurate  temperature  prediction  at 
35  mm,  but  the  prediction  of  OH  is  not  satisfactory  and  there  is 
a  large  discrepancy  between  OH  modeling  results  and 
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Fig.  8.  Effect  of  turbulent  combustion  models  (adopting  the  modified  SKE  turbulence  model,  GRI-Mech  3.0  mechanism  and  (DO)  radiation  model):  comparison  between  exper¬ 
imental  measurement  [16]  and  model  predictions  (temperature  (a,  c,  e)  and  OH  number  density  (b,  d,  f))  for  3%  02  and  9%  02  ethylene/nitrogen  flames. 


measurements.  The  modeling  results  indicate  that  the  flame  is  lif¬ 
ted  and  that  there  are  no  reactions  at  the  location.  Experimentally, 
however,  despite  the  visual  appearance  of  a  lifted  flame,  a  reaction 
zone  was  indeed  recorded  in  this  region.  For  the  3%  O2  flame  at 
125  mm,  the  models  generate  a  reasonable  prediction  of  temper¬ 
ature  and  OH  number  density.  The  9%  O2  in  the  coflow  stream 
makes  this  flame  more  reactive  than  the  equivalent  3%  case. 
Nonetheless,  these  flames  still  visually  appeared  lifted  during  the 
experiments.  This  effect  is  largely  overestimated  by  EDC,  modified 
EDC  and  flamelets  models,  while  it  is  correctly  predicted  by  the  PDF 
transport  models. 

Overall,  the  results  demonstrate  that  the  modified  EDC  and  PDF 
transport  (EMST  and  MC)  models  are  adequate  for  temperature 
prediction  in  the  3%  O2  flames  of  all  four  fuels  studied.  For  the  3%  O2 
flames,  all  models  under-predict  the  magnitude  of  the  OH  peak,  but 
the  PDF  transport  models  give  a  better  prediction  of  location  of  OH 
concentration  peak.  In  addition,  the  values  of  OH  number  density  in 
the  coflow  stream  predicted  by  PDF  transport  models  are  more 
satisfactory  than  the  other  models  for  3%  flames.  For  the  9%  O2 


flame  of  different  fuels,  the  PDF  transport  models  yield  reasonably 
good  agreement  with  experiments  for  the  temperature  and  OH 
number  density. 

4.2.5.  Effect  of  fuel  composition 

Fig.  9  shows  a  comparison  between  the  measured  peak  values  of 
temperature  and  OH  concentration  and  the  corresponding  pre¬ 
dictions  of  the  models  for  the  different  flames.  In  this  way,  this 
figure  summarizes  the  results  already  presented  in  more  detail  in 
the  previous  figures,  but  allows  the  trends  associated  with  the 
changes  in  fuel  composition  to  be  highlighted.  It  is  interesting  to 
observe  that  at  a  distance  of  35  mm  from  the  inlet,  the  maximum 
temperature  is  a  good  indicator  of  the  reactivity  of  the  system.  The 
experimental  values  reported  in  Fig.  9a  clearly  show  that  flame 
number  3  (partially  premixed  C2H4/air)  is  the  most  reactive.  All  the 
models  are  able  to  reproduce  this  trend.  The  PDF  transport  models 
significantly  underestimate  the  temperature  peak  for  flame  3, 
whilst  giving  satisfactory  predictions  for  the  other  flames.  It  is  also 
possible  to  observe  that  the  EDC  model  overestimates  the  reactivity 
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Fig.  9.  Maximum  temperature  and  OH  number  density  of  experimental  measurement  [16]  and  model  predictions  (different  turbulent  combustion  models)  for  studied  flames.  Flame 
number  1:  (C2H4),  flame  number  2:  (C2H4/H2),  flame  number  3:  (C2H4/air)  and  flame  number  4:  (C2H4/N2). 


for  flame  1  (C2H4)  and  flame  2  (C2H4/H2),  while  it  is  in  better 
agreement  with  measured  values  for  the  diluted  C2H4/N2  flame.  As 
already  observed,  the  modified  EDC  provides  better  predictions  in 
these  conditions.  This  is  not  surprising  since  this  modification  of 
the  EDC  model  was  proposed  by  De  et  al.  [15]  for  methane  under 
MILD  combustion  conditions.  In  this  work  we  extended  the  vali¬ 
dation  of  the  model  to  ethylene  mixtures.  The  flamelets  model  does 
not  satisfactorily  predict  the  trend  of  the  temperature  peak  with 
the  fuel  composition,  especially  for  the  ethylene/air  and  ethylene/ 
hydrogen  flames.  As  already  discussed,  the  OH  concentration  is 
under-estimated  by  the  models,  which  are  not  able  to  explain 
the  large  OH  peak  of  the  ethylene/hydrogen  flame  at  35  mm.  The 
results  at  125  mm  could  be  affected  by  the  entrainment  of  the 
surrounding  air  into  the  flame,  which  is  difficult  to  accurately 
capture  in  the  models.  Some  of  the  deviations  presented  in  Fig.  9 
can  be,  at  least  partially,  explained  on  this  basis.  Fig.  9  also  shows 
the  comparison  of  temperature  peaks  and  OH  concentrations  for 
the  9%  O2  flames  (Fig.  9e  and  f).  It  is  evident  that  only  the  PDF 
transport  models  (MC  and  EMST)  are  able  to  characterize  the 
reactivity  of  the  system  in  these  conditions.  The  ethylene/hydrogen 


flame  shows  the  highest  temperature  peak,  while  the  ethylene/air 
flame  shows  a  relatively  low  temperature  peak  compared  with  the 
3%  O2  case.  The  EDC  and  the  flamelet  model  do  not  correctly 
describe  the  system’s  reactivity  and  predict  a  shorter  liftoff  height, 
therefore  over-estimating  the  temperature  peak.  The  same  analysis 
applies  to  the  OH  peaks. 

To  better  analyze  the  effect  of  the  turbulent  combustion  model 
on  the  OH  predictions,  Fig.  10  shows  the  OH  distribution  in  the 
flame  region  close  to  the  burner  exit  plane.  This  comparison 
confirms  that  the  EDC  model  predicts  an  attached  flame  for  all  the 
conditions,  while  the  PDF  model  shows  a  significant  liftoff  height 
for  the  ethylene  flames  diluted  with  air  or  nitrogen.  This  different 
behavior  is  clearly  evident  in  the  OH  predictions  discussed  in  Fig.  9. 
It  is  interesting  to  note  that  experimentally  the  apparent  liftoff 
for  the  9%  02  flames  is  26  mm  for  the  C2H4  flame,  33  mm  for  the 
C2H4/Air  flame  and  34  mm  for  the  C2H4/N2  flame.  This  behavior 
is  quite  well  reproduced  by  the  PDF  model.  On  the  contrary,  the 
C2H4/H2  flame  is  attached  and  this  is  not  completely  captured  by 
the  model.  The  3%  02  flames  are  all  weakly  reactive  but  all  attached 
to  the  burner,  as  experimentally  observed.  Both  the  EDC  and  PDF 
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Fig.  10.  Predicted  OH  number  densities  for  the  studied  flames  using  the  EDC  and  PDF-EMST  models. 


models  correctly  predict  this  behavior,  based  on  the  analysis  of  the 
OH  plots.  Medwell  et  al.  [16]  showed  that  rather  than  referring  to 
them  as  lifted,  these  ethylene  flames  can  be  better  described  as 
having  a  transition  in  the  reaction.  The  apparent  liftoff  height  cor¬ 
responds  to  a  transition  of  weak  to  strong  OH  levels.  As  such,  rather 
than  identifying  a  “liftoff’  height,  it  is  more  appropriate  to  refer 
to  this  as  a  transition  point.  The  general  observations  relating  to 
this  transition  point  seem  to  be  analogous  to  conventional  lifted 
flame  zones. 

Table  2  reports  the  estimated  error  of  the  predicted  apparent 
liftoff  heights  (transition  points),  for  the  different  turbulent  com¬ 
bustion  models.  This  deviation  is  simply  calculated  as  the  difference 
(in  axial  space)  between  the  35  mm  location  and  the  region  where 
the  radial  predicted  temperature  profile  matches  the  one  measured 
at  35  mm.  As  an  example,  for  the  ethylene  flame,  the  EDC  model 
largely  over-predicts  the  temperature  profile  at  35  mm,  especially 
for  the  9%  02  flame  (see  Fig.  5a  and  e).  For  this  reason,  to  match  the 
measured  and  predicted  temperature  profiles,  it  is  necessary  to 
shift  the  axial  position  toward  the  inlet  by  25  mm  ( i.e .  using  the 


Table  2 

Estimated  apparent  liftoff  (transition  points)  error  (ELOE)  (units:  mm)  of  different 
turbulent  combustion  models.  The  estimate  is  based  on  the  comparison  of  predicted 
and  measured  temperatures. 


c2h4 

c2h4/h2 

C2H4/air 

c2h4/n2 

02  in  coflow 
Experiment 

3% 

9% 

3% 

9% 

3% 

9% 

3%  9% 

EDC 

-25 

-30.8 

-27.5 

-31.1 

-9 

-30.7 

-9  -28.6 

Modified  EDC 

-17 

-29 

-16 

-23 

+6 

-29.9 

-8  -27.6 

Flamelet 

-18 

-28.5 

-21 

-19 

+9 

-33 

-5  -30.8 

PDF  transport 

-8 

-19.5 

-4 

-12 

+6 

+4 

+3  -1 

(Modified  curl) 

PDF  transport  (EMST) 

-7 

-18.5 

-5 

-15 

+5 

0 

-2  -0.9 

profile  at  the  10  mm  axial  location)  for  the  3%  case  and  of  30.8  mm 
for  the  9%  case.  The  shift  is  larger  for  the  9%  case  because  the 
temperature  is  largely  overestimated  by  the  model.  We  use  this 
shift  as  an  indicator  of  the  error  in  the  transition  points  of  the 
flames.  It  is  important  to  notice  that  this  indicator  shows  a  devia¬ 
tion  which  can  be  due  both  to  an  error  in  the  predicted  liftoff  height 
and  in  the  reactivity  of  the  flame.  Since  they  are  all  attached,  the 
negative  values  of  ELOE  for  the  3%  O2  flames  are  likely  to  indicate 
that  the  EDC  model  in  particular  overestimates  the  reactivity  of  the 
flame  and  the  position  of  the  transition  point  between  weak  and 
strong  reactivity.  Also  in  these  cases,  the  PDF  models  are  in  better 
agreement  with  the  experiments. 

The  comparison  shows  the  better  agreement  of  the  PDF  trans¬ 
port  models  compared  with  EDC  and  flamelet  models,  which  is 
likely  due  to  their  ability  to  characterize  localized  extinction  or 
intermittency  phenomena.  Another  observation  is  that  the  modi¬ 
fied  EDC,  compared  to  the  original  model,  improves  the  predictions 
which  become  very  similar  to  the  ones  of  the  flamelet  model. 
Finally,  the  errors  are  larger  for  the  9%  02  cases  relative  to  3%  O2. 

4.3.  Formation  and  consumption  ofCH20 

Formaldehyde  is  an  important  intermediate  in  the  oxidation  of 
hydrocarbons.  It  is  also  formed  in  significant  amounts  in  the  low 
temperature  combustion  of  hydrocarbons  and  oxygenated  biofuels. 
It  is  interesting  to  note  that  the  premixing  of  O2  with  the  fuel  allows 
CH20  to  be  formed  in  regions  of  the  flame  where  the  local  condi¬ 
tions  make  its  consumption  a  slow  process.  The  expected  effect 
of  premixing  is  to  increase  the  peak  of  formaldehyde  in  the  flame, 
up  to  five  times  as  compared  to  purely  diffusive  flames  [46].  A 
detailed  analysis  of  CH20  formation  in  flames  in  a  heated  and 
diluted  oxidant  stream  has  been  presented  by  Medwell  et  al.  [47] 
using  laminar  flame  calculations.  They  showed  that  under  MILD 


550 


S.R.  Shabanian  et  al.  /  Applied  Thermal  Engineering  52  (2013)  538-554 


conditions  O2  penetrates  the  reaction  zone,  leading  to  partial  pre¬ 
mixing  and  thus  contributing  to  the  improved  flame  stability  under 
these  conditions. 

Formaldehyde  profiles  were  measured  by  Medwell  et  al.  [16]  and 
provided  in  arbitrary  units  for  the  JHC  burner.  Since  the  experimental 
measurements  are  not  quantitative,  we  decided  to  apply  a  scaling 
factor  to  the  predicted  CH2O  number  density.  For  each  condition 
(figure),  the  scaling  factor  is  calculated  to  match  the  measured  and 
predicted  CH2O  peak  in  the  pure  ethylene  flame.  The  same  factor  is 
then  applied  to  the  predicted  profiles  of  the  other  flames.  In  this  way, 
it  is  possible  to  compare  the  shape  of  the  profiles  and  the  effect  of 
the  fuel  on  the  formation  of  formaldehyde.  Fig.  11  shows  that  the 
modified  EDC  and  the  PDF  models  are  able  to  correctly  characterize 
the  relative  formation  of  CH2O  in  the  pure  ethylene  and  in  the  par¬ 
tially  premixed  flames.  The  relatively  large  amount  of  formaldehyde 


on  the  axis  of  the  flame  is  not  predicted  by  the  models.  Fig.  Ilf 
summarizes  the  results  in  terms  of  peak  CH2O.  Due  to  the  normal¬ 
ization,  all  the  predicted  CH20  peaks  coincide  with  the  measured 
peak  CH2O  value  for  the  C2H4  flame.  It  is  quite  evident  that  all  models 
succeeded  in  predicting  the  effect  of  fuel  composition  on  the  for¬ 
mation  of  CH2O.  It  is  interesting  to  observe  that  the  effect  of  H2  is  to 
reduce  the  formation  of  CH2O.  This  is  expected  as  the  availability  of  C 
atoms  is  reduced  due  to  the  dilution  with  hydrogen.  The  same  effect 
can  be  observed  for  N2  dilution.  On  the  contrary,  the  partial  pre¬ 
mixing  with  air  (O2)  has  a  relevant  chemical  effect  which  is  properly 
captured  by  the  models.  To  compare  the  effect  of  the  kinetic  mech¬ 
anism,  Fig.  12  shows  the  predicted  CH2O  formed  in  two  counterflow 
laminar  flames  calculated  in  conditions  representative  of  those  of  the 
experiments.  The  oxidant  has  the  temperature  and  composition  of 
the  coflow  used  in  the  JHC  flames  and  the  fuel  is  pure  ethylene  or 
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Fig.  11.  Experimental  [16]  and  modeling  mean  radial  profile  of  formaldehyde  for  3%  flame  at  axial  location  of  35  mm.  Orange:  C2H4,  blue:  C2H4/H2,  black:  C2H4/air  and  green:  C2H4/ 
N2.  Lines:  experimental  results  and  dashed  lines:  modeling  results,  (panels  a-e).  Panel  f:  comparison  between  measured  and  predicted  maximum  CH20  (of  panels  a-e)  using 
different  turbulent  combustion  models.  Flame  number  1 :  (C2H4),  flame  number  2:  (C2H4/H2),  flame  number  3:  (C2H4/air)  and  flame  number  4:  (C2H4/N2).  (For  interpretation  of  the 
references  to  colour  in  this  figure  legend,  the  reader  is  referred  to  the  web  version  of  this  article.) 
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Fig.  12.  Comparison  between  formaldehyde  mole  fraction  of  two  kinetic  mechanisms  (GRI-Mech  3.0  and  POLIM1)  found  from  laminar  flame  calculation.  9%  02  coflow  oxidant 
stream.  Tox|  =  1100  K.  a  ~  200  s 


ethylene/air.  It  is  evident  that  the  formation  of  CH2O  in  the  partially 
premixed  case  is  significantly  larger  and  that  it  occurs  earlier  (closer 
to  the  fuel  inlet),  while  in  the  diffusive  flame  CH2O  is  mainly  formed 
in  the  lean  side  of  the  flame  front.  The  CH20  predicted  by  the 
GRI-Mech  3.0  model  is  larger  than  the  amount  predicted  by  the 
POL1M1  model  in  the  partially  premixed  flame.  The  reaction  flux 
analyses  presented  in  Figs.  13  and  14  clearly  show  the  difference 
between  ethylene  premixed  and  diffusion  flames.  Both  for  POLIMI 
and  GRI-Mech  3.0,  the  consumption  of  ethylene  proceeds  through 
the  formation  of  vinyl  radicals  which  subsequently  form  acetylene. 
The  fate  of  vinyl  radicals  is  different  when  O2  is  available  in  the 
system  at  the  early  stage  of  ethylene  consumption.  Vinyl  radicals 
react  with  O2  to  produce  CH2CHO,  which  subsequently  forms  CH2CO. 
Ketene  then  forms  formaldehyde  either  directly  (POLIMI)  or  via  CH3 
formation  and  oxidation.  The  comparison  of  Figs.  13  and  14  also 
explains  the  differences  between  partially  premixed  and  diffusive 
flames.  It  is  even  more  pronounced  at  the  125  mm  axial  location 
(compared  with  the  35  mm  location)  that  the  partially  premixed 
flame  is  characterized  by  the  largest  amount  of  CH2O  (refer  to  Fig.  S4 
(Supplementary  data)  for  further  details).  The  effect  of  fuel  compo¬ 
sition  is  quite  well  captured  by  the  combustion  models.  Only  the 
flamelet  model  failed  in  predicting  this  observed  value  for  the  par¬ 
tially  premixed  flame.  For  the  9%  O2  flame  the  previous  findings  are 


confirmed  at  35  mm:  the  PDF  transport  and  modified  EDC  models 
provide  a  satisfactory  agreement  with  the  experimental  data  while 
the  flamelets  model  shows  a  large  deviation  especially  in  the  parti¬ 
ally  premixed  conditions  (refer  to  Fig.  S5  (Supplementary  data)  for 
further  details).  Also  in  this  case  the  modified  EDC  improves  the 
agreement  with  the  measurements.  In  fact,  the  EDC  overestimates 
the  increase  of  CIT2O  when  air  is  premixed  with  the  fuel.  This  devi¬ 
ation  was  not  present  in  the  3%  O2  flame. 

4.4.  Effect  of  the  kinetic  mechanism 

The  effect  of  the  kinetic  mechanism  on  the  modeling  was 
investigated  by  comparing  predictions  using  two  different  mech¬ 
anisms,  POLIMI  and  GRI-Mech  3.0.  First  we  discuss  the  difference 
between  the  two  mechanisms  in  laminar  flame  conditions,  then  in 
the  ethylene  turbulent  flames  studied  in  this  work. 

Fig.  15a  shows  a  comparison  between  experimental  measure¬ 
ments  [48—52]  and  model  predications  of  laminar  flame  speed  of 
ethylene  in  air  for  the  kinetic  mechanisms  considered  in  this  study. 
The  flame  speeds  were  calculated  using  the  OpenSMOKE  code  [53]. 
A  more  detailed  description  on  the  flame  speed  calculations  can  be 
found  in  Ref.  [34],  which  also  contains  an  extensive  validation  of 
the  POLIMI  model.  Two  observations  can  be  made:  first,  it  is  quite 
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Fig.  13.  a.  Reaction  path  analysis  for  GRI-Mech  3.0  and  the  conditions  of  Fig.  12a  (diffusion  flame),  b.  Reaction  path  analysis  for  GRI-Mech  3.0  and  the  conditions  of  Fig.  12b  (partially 
premixed  flame). 
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Fig.  14.  a.  Reaction  path  analysis  for  POL1M1  and  the  conditions  of  Fig.  12a  (diffusion  flame),  b.  Reaction  path  analysis  for  POL1M1  and  the  conditions  of  Fig.  12b  (partially  premixed  flame). 


evident  the  larger  reactivity  of  ethylene  compared  to  methane, 
whose  maximum  flame  speed  in  these  conditions  is  35—38  cm/s 
under  these  conditions.  Second,  it  can  be  observed  that  the 
reduced  C2H4  model  developed  at  Politecnico  di  Milano  is  able  to 


equivalence  ratio 


Fig.  15.  a.  Laminar  flame  speed  of  ethylene  in  air  at  1  atm  and  300  K.  Comparison 
between  GRl-Mech  3.0  and  POL1M1  predicted  values  with  experimental  data  [48-52],  b. 
Laminar  flame  speed  of  ethylene  in  hot  diluted  air  (9%  02)  at  1  atm  and  initial  tem¬ 
perature  of  1000  K.  Comparison  between  GRl-Mech  3.0  and  POL1M1  predicted  values. 


correctly  predict  the  reactivity  of  ethylene  and  the  velocity  of  flame 
propagation.  On  the  other  hand,  GRl-Mech  3.0  —  a  mechanism 
which  was  originally  developed  for  methane  and  natural  gas 
combustion  —  significantly  overestimates  the  flame  speed.  The 
same  difference  in  the  reactivity  of  the  flame  can  be  observed  in 
Fig.  15b,  where  the  flame  speed  is  calculated  for  a  mixture  of 
ethylene  and  hot  and  diluted  air.  The  oxidizer  contains  9%  O2  and 
has  the  same  composition  of  the  oxidizer  used  in  the  coflow  of  the 
JHC  burner  experiments  discussed  in  this  paper.  Referring  back  to 
the  earlier  discussion  surrounding  liftoff  height,  it  may  be  deduced 
that  the  liftoff  behavior  was  not  captured  in  the  model  also  as 
a  result  of  the  excessively  high  flame  speed  when  using  the 
GRl-Mech  3.0  mechanism. 

The  EDC  combustion  model  and  modified  SKE  turbulence 
model  have  been  adopted  in  simulations  of  the  9%  O2  ethylene 
flames  and  show  that  the  temperature  and  OH  number  density 
profiles  are  very  similar  for  these  two  mechanisms  (refer  to 
Fig.  S6  (Supplementary  data)  for  further  details).  It  is  interesting 
to  note  that,  despite  differences  in  the  laminar  flame  speed  for  the 
two  chemistry  models  considered,  the  choice  of  kinetics  does  not 
significantly  influence  the  reaction  zone  profiles  in  the  turbulent 
flame  conditions.  This  suggests  that  the  turbulence— chemistry 
interactions  play  a  major  role  that  largely  supersedes  the  differ¬ 
ences  in  the  kinetic  schemes.  On  the  contrary,  Fig.  16  shows  that 
the  tendency  of  GRl-Mech  3.0  to  form  a  larger  amount  of  CH2O  is 
also  evident  in  the  turbulent  partially  premixed  flame.  For  the  3% 
O2  coflow  cases,  it  is  difficult  to  stabilize  the  flame  with  the 
POLIM1  mechanism  and  the  EDC  turbulent  combustion  model. 
This  effect,  which  prevents  a  wide  range  of  comparisons  using  the 
POLIM1  kinetic  scheme,  may  be  due  to  the  already  discussed 
lower  reactivity  of  the  POLIMI  mechanism  with  respect  to  the 
GRl-Mech  3.0,  and  to  the  difficulty  of  EDC  in  describing  the 
mechanism  of  jet  flame  stabilization,  as  already  noted  by 
Zahirovic  et  al.  [54], 

It  is  important  to  note  that  although  PDF  transport  models 
provide  simulation  results  which  are  in  better  agreement  with 
experimental  measurements  the  EDC  remains  a  valuable  model. 
The  interest  in  EDC  is  mainly  motivated  by  the  advantages  related 
to  its  relatively  low  computational  requirements.  This  is  important, 
especially  when  referring  to  CFD  modeling  of  real  three- 
dimensional  combustors.  The  complexity  and  the  size  of  the  com¬ 
putational  grid  associated  with  such  a  geometry  makes  it  difficult 
to  use  PDF  transport  models  with  detailed  chemistry. 
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Fig.  16.  Experimental  and  modeling  mean  radial  profile  of  formaldehyde  for  9%  flame  at  axial  location  of  35  mm.  Orange:  C2H4,  blue:  C2H4/H2,  black:  C2H4/air  and  green:  C2H4/N2. 
Lines:  experimental  results  and  dashed  lines:  modeling  results  obtained  using  the  EDC  model  and  two  different  kinetic  mechanisms.  (Panels  a  and  b).  Predicted  and  measured 
maximum  CH20  using  two  kinetic  mechanisms.  Flame  number  1:  (C2FL4),  flame  number  2:  (C2H.4/H2),  flame  number  3:  (C2H4/air)  and  flame  number  4:  (C2H4/N2)  (panel  c).  (For 
interpretation  of  the  references  to  colour  in  this  figure  legend,  the  reader  is  referred  to  the  web  version  of  this  article.) 


5.  Conclusions 

Turbulent  nonpremixed  jet  flames  issuing  into  a  heated  and 
highly  diluted  coflow  were  investigated  in  this  study  using  a  CFD 
code  (FLUENT  12.1.4).  An  extensive  set  of  simulation  results  is  pre¬ 
sented  for  several  ethylene  flames,  either  undiluted  or  diluted  with 
hydrogen,  air  or  nitrogen,  to  highlight  the  effect  of  using  different 
turbulent  combustion  models  and  kinetic  schemes  (GRi-Mech  3.0 
and  POLIMI).  The  predictions  of  the  flame  structure  using  the  two 
kinetic  mechanisms  are  very  similar.  In  the  case  of  the  3%  O2  flames, 
using  the  POLIMI  mechanism  with  the  EDC  model  led  to  flame 
extinction.  This  behavior  is  attributed  to  the  difficulties  of  the  EDC 
model  in  predicting  the  jet  flame  stabilization  [54],  The  modified  k—e 
model  succeeded  in  simulating  the  fuel  jet  fluid  dynamics  flame 
structure.  Thermal  radiation  was  included  in  the  model  using  the 
discrete  ordinate  model.  The  comparison  between  the  predicted 
and  the  measured  mean  temperature  profiles  showed  a  systematic 
deviation  for  the  standard  EDC  model.  This  model  is  found  to  over¬ 
estimate  the  reactivity  of  the  flame,  both  in  the  3%  and  9%  02  cases. 
In  other  words,  the  EDC  tends  to  under-predict  liftoff  height.  The 
predictions  are  significantly  improved  by  using  a  modified  version  of 
the  EDC  [15,24],  which  was  developed  for  MILD  combustion  condi¬ 
tions.  The  modeling  results  demonstrate  that  the  modified  EDCyields 
reasonable  results  with  a  relatively  low  computational  effort.  The 
agreement  between  measurements  and  predictions  of  the  PDF 
transport  and  modified  EDC  models  is  generally  better  than  for  the 
flamelet  model.  A  comparative  study  is  also  presented  for  the  impact 
of  two  micro-mixing  models  in  transported  scalar  PDF  simulations  of 
these  flames.  The  micro-mixing  models  are  MC  (modified  curl)  and 
EMST  (Euclidean  minimum  spanning  tree).  In  general  the  difference 


between  the  predictions  of  the  two  PDF  models  is  not  significant. 
The  predictions  are  in  good  agreement  with  the  measurements,  and 
are  improved  especially  for  the  apparent  liftoff  heights  and  the 
temperature  peaks.  The  effect  of  the  fuel  composition  on  the  CH2O 
formation  is  also  discussed  and  shows  good  agreement. 
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